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Exciton effects on soliton and bipolaron lattice states are investigated us- 



ing an electron-lattice Peierls model with long-range Coulomb interactions. 
The Hartree-Fock (HF) approximation and the single-excitation configuration- 
interaction (single-CI) method are used to obtain optical absorption spectra. 
We discuss the following properties: (1) The attraction between the excited 
electron and the remaining hole makes the excitation energy smaller when the 
correlations are taken into account by the single-CI. The oscillator strengths of 
the lower excited states become relatively larger than in the HF calculations. 

(2) We look at variations of relative oscillator strengths of two or three kinds of 
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excitons described by the single-CI. While the excess-electron concentration is 
small, the ratio of the oscillator strengths of the exciton with the lowest energy, 
which is calculated against the total electronic excitation oscillator strengths, 
increases almost linearly. The oscillator strengths accumulate at this exciton 
as the concentration increases. 



PACS numbers: 7135, 7840, 7138 
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1. Introduction 



It has been well known that correlation effects are present among 7r-electrons 
in conjugated polymers. For example: 

• The envelope of the wavefunction of the midgap state of the neutral 
(spin) soliton in trans-polyacetylene has sites with spin density of the 
opposite sign [1]. This fact of the "negative spin density" has been ex- 
plained by considering Coulomb interactions in the Su-Schrieffer-Heeger 
(SSH) model [2] of conjugated polymers [3]. 

• The electron paramagnetic resonance experiment of pernigraniline base 
[4], reported recently, has studied the spin distribution around a neutral 
soliton. It has been suggested that sites with negative spin density exist 
as the consequence of correlation effects. 

• Nonlinear optical response functions of polydiacetylenes exhibit excita- 
tion structures owing to the presence of excitons. Their structures have 
been explained theoretically by using the intermediate exciton formalism 
[5]. 

Therefore, it is generally interesting to study Coulomb interaction effects in 
conjugated polymers. In fact, this problem has been considered by various 
authors for more than a decade since the SSH model [2] was proposed. 

Electronic excitation structures in the half-filled conjugated polymers with 
the constant dimerization have been theoretically investigated by using the 
exciton formalism [5] and the exact diagonalization method [6], and also by 
solving the time- dependent Hartree-Fock (HF) equations [7]. This problem 
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was pointed out earlier, but excitation structures have been considered inten- 
sively only recently, relating with origins of nonlinear optical spectra [6-8]. 
The most remarkable consequence of correlation effects is that the lowest en- 
ergy excitation has the largest oscillator strength. It is clearly seen when 
the optical spectra calculated by using the HF wavefunctions are compared 
with the spectra with the correlation effects. This fact does not depend on 
whether the higher correlations are taken into account by the single-excitation 
configuration-interaction (single-CI) method [5] , or by the time-dependent HF 
formalism [7]. 

It is widely known that the soliton, polaron, and bipolaron lattices are 
present [9], when the SSH model [2], its continuum version, and the extended 
model with the term of the nondegeneracy [10] are doped with electrons or 
holes. New bands related with nonlinear excitations develop in the Peierls gap 
as the doping proceeds. When correlation effects are considered by the single- 
CI, the excitation structures exhibit the presence of excitons. There is one 
kind of exciton in the half-filled system, where the excited electron (hole) sits 
at the bottom of the conduction band (top of the valence band). We will call 
this as the "intercontinuum exciton" . In the soliton lattice states of the doped 
SSH model for degenerate conjugated polymers, there are small gaps between 
the soliton band and the continuum states, i.e., valence and conduction bands. 
Therefore, the number of the kind of excitons would increase and their presence 
will be reflected in structures of the optical spectra. A new exciton, namely, 
"soliton-continuum exciton" will appear when the electron-hole excitation is 
considered between the soliton and one of the continuum bands. The main 
purpose of this paper is to embody the above picture. There have been a lot 
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of investigations of correlation effects in doped SSH systems, but the work in 
the view point of excitons has been rarely performed. General properties of 
exciton effects on soliton lattice systems are the first interest of the paper. The 
portion of the spectral weights of the two kinds of excitons will be calculated 
and discussed. 

There are two midgap bands, when the nondegenerate conjugated polymers 
are doped and bipolaron lattice states are formed. The excitation spectra will 
become complicated due to the increase of number of kinds of excitons. The 
second part of this paper will be devoted to this problem, and the relative 
spectral weight of each exciton will be studied again. 

This paper is composed as follows. In §2, the model is introduced and 
the numerical method is explained. Results for the soliton lattice system are 
reported in §3. Effects of the nondegeneracy are investigated in §4. The paper 
is summarized in §5, and discussion is given the final section. 



The following hamiltonian is used to discuss excitonic effects in soliton and 
bipolaron lattice states of Peierls systems: 



2. 



Model 



H = H sm + H mt . 



(1) 



The first term -£/ssh of Eq. (1) is the SSH model: 



#ssh = - ®Vi + (- 1 )%t\(cl a c i+1>IT + h.c.) 



(2) 
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where t is the hopping integral of the system without the dimerization; a is the 
electron-phonon coupling constant which changes the hopping integral linearly 
with respect to the bond variable yf, Sot is the Brazovskii-Kirova (BK) term 
which measures the degree of the nondegeneracy (it was originally introduced 
in the continuum model [10]); Cj jCr is an annihilation operator of the 7r-electron 
at the site % with spin a; the sum is taken over all the lattice sites of the 
periodic chain; and the last term with the spring constant K is the harmonic 
energy of the classical spring simulating the a-bond effects. The second term 
of Eq. (1) is the long-range Coulomb interaction in the form of the Ohno 
potential [11]: 

= U XX C !,T C *,T ~~ -^-)( c li c hl ~ -f-) ( 3 ) 
+ 52 w ( r i,j) E c l,aCi,a - n el ) (52 c\ T c j}T - n e i) , 

where n e \ is the number of 7r-electrons per site, r it j is the distance between the 
ith and jth sites, and 

W(r) = . 1 (4) 

^/(i/uy + (r/r vy 

is the Ohno potential. The quantity W(0) = U is the strength of the onsite 
interaction, and V means the strength of the long range part. 

The model is treated by the HF approximation and the single-CI for the 
Coulomb potential. The adiabatic approximation is forced on the bond vari- 
ables. The HF order parameters and bond variables are determined selfconsis- 
tently using the standard iteration method [12]. After the HF approximation 
H =^ H UF , we divide the total hamiltonian as H = H UF + H'. The term H' 
becomes 

H' = U 52(4,1*,! -A,t)(4,| c U -Pi,l) ( 5 ) 



+ ]T W(r itj )[J2(cl a Ci, a - p it a)(cl T C jtT - p jtT ) 

+ 5Z( r i,j,crCj/ i0 -Cj i0 - + T j,i,CT C l t a C j,a ~ T i,j,v T j,i,a)] i 
a 

where p ij(T = (cj )0 .Cj )0 .) and t^^ c = (c\ a Cj ja ) are Hartree-Fock order parameters. 
When we write the Hartree-Fock ground state \g) = IlA:occupicd c a,t c a,J0) an d 
the single electron-hole excitations \p\) = c^ a c\ jT \g) (/i means an unoccupied 
state; we assume both singlet and triplet excitations in this abbreviated no- 
tation), the matrix elements of the HF part and the excitation hamiltonian 
become as follows: 

(//A'|(#hf-(#hf))|M) = V,^v,a(£ m -£a), (6) 
(p'X'\(H'-{H'))\p\) = 2JS S -K, (7) 

where (. . .) means the expectation value with respect to the HF ground state, 
i.e., (H HF ) = (g\H HF \g) and (H') = (g\H'\g), and is the energy of the HF 
orbital, 5s = 1 for spin singlet, 5s = for spin triplet, and 

J(//,A';//,A) = J2 V ijW)W)mm (8) 
K(p',\';p,\) = J2 V i,Mj)W)(M(i\V (9) 

with V i:i = U, Vij = W(r it j) for i ^ j. The diagonalization of the total 
hamiltonian H gives the set of the excited states {\k)} within the single-CI 
method. In the actual calculation, we limit the spin configurations to the 
singlet excitations which are the main interests of optical excitations. 

We assume a geometry of a ring for a polymer chain, in order to remove 
edge effects. If we use an open boundary, the dimerization and thus the Peierls 
gap becomes larger near the two edges, and this might result in artifacts of 
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optical spectra. We shall use the coordinate of jth carbon atoms, 
2nj . 2nj 

(rcos— ,rsm— ,0), (10) 

where r = Na/ (2n) is the radius of the polymer ring; N is the system size 
and a is the lattice constant. The electric field of light is parallel to the x-y 
plane. In order to obtain optical spectra which are independent of the relative 
positions of solitons with respect to the direction of light, we shall sum up 
two spectra where light is along with the x- and y-directions. So, we use the 
following formula of the spectrum: 

Y,E k P{uj - E K )((g\x\K}(K\x\g) + (g\y\ K ) (n\y\g)). (11) 

K 

Here, P(lo) = ^/[n(uj 2 + r ) 2 )} is the Lorentzian distribution (7 is the width), 
E K is the electron-hole excitation energy, and \g) means the ground state. In 
eq. (11), the quantity, 

E K ((g\x\K)(K\x\g} + (g\y\K)(K\y\g)), (12) 

is the oscillator strength of the excited state \k). Applying electric field to the 
ring-shaped polymer simulates the situation that polymer chains are oriented 
randomly in every direction within the x-y plane. The average over orienta- 
tions are effectively performed. The similar idea has been used in the recent 
paper by Abe et al [5] . 

The system size is chosen as N = 80, 100, 120 when the electron number 
is even (it is varied from N el = N, N + 2, N + 4, N + Q to iV + 8), because 
the size around 100 is known to give well the energy gap value of the infinite 
chain. More larger system size becomes tedious for doing single-CI calcula- 
tions which call for huge computer memories. When there is one soliton, we 



take N = 81, 101, 121, and use the periodic boundary condition also. Both the 
combinations, (N,N C \) = (80,84) and (120, 126), give the 5% concentration of 
the excess electrons. Therefore, two symbols will appear at the 5% concen- 
tration in Figs. 4, 5, 7, and 8, when numerical results are plotted against the 
doping concentration. This is a remark for later convenience. 

In principle, we have to adjust parameters and find appropriate ones in 
order to reproduce experimental data, such as, the energy gap and the dimer- 
ization amplitude. But, we will change parameters arbitrary in a reasonable 
range in order to look at excitonic effects clearly. The Coulomb parameters 
are changed within < V < U < 5t, and we show results for U = 2V = 2t and 
= At as the representative cases. Other parameters, t = 1.8eV, K = 21eV/A 2 , 
and a = 4.1eV/A, are fixed in view of the general interests of this paper. All 
the quantities of energy dimension are shown in the units of t. 

3. Soliton lattice systems 

Figure 1 shows the typical lattice configuration and excess-electron density 
distribution for N = 100, iV el = 104, U = 4t, V = 2t, and 5 = 0. Both 
quantities are the smoothed data after the removal of small oscillations between 
even- and odd-number sites. There are four charged solitons due to the excess- 
electron number N c \ — N = 4. Solitons are arrayed equidistantly. The excess- 
electron density shows the oscillation of the charge density with its maxima 
at the soliton centers. 

Next, let us look at optical spectra and consider exciton effects. We calcu- 
late also the optical spectra from single-electron excitations among HF energy 
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levels as references of exciton effects. The typical optical spectra within the 
HF approximation and with the single-CI are shown as Figs. 2 and 3, respec- 
tively. The broadening 7 = 0.05t is used. The Coulomb parameters are U = At 
and V = 2t. Relatively strong Coulomb interactions are taken in order to look 
at the exciton effects in the optical response clearly. The system size and the 
electron number are (N,N el ) = (101, 102), (100, 102), (100, 104), for (a), (b), 
and (c) (of Figs. 2 and 3), respectively. 

In order to characterize properties of optical excitations, in other words, 
to identify whether the excitons are formed between valence and conduction 
bands ( "intercontinuum exciton"), or between soliton and conduction bands 
( "soliton-continuum exciton"), we calculate the component of the electron- 
hole excitation between the soliton and conduction bands for the absorption 
spectra within HF as well as those with the single-CI. Here, in other words, the 
"component" is the ratio of the oscillator strengths of the soliton-continuum 
exciton with respect to the total oscillator strengths of the two kinds of exci- 
tons. Results are shown in Figs. 2 and 3, superposed with the optical spectra. 
We can easily determine the position of the optical gap by comparing compo- 
nents of each kind of exciton. The optical gap of one exciton is located at the 
lowest energy where the component is larger than that of the other exciton. 

As we proceed from Fig. 2(a) to Fig. 2(c) with increasing the soliton 
concentration, the contribution from the transition between the soliton band 
and the conduction band becomes larger than that between the valence and 
conduction bands. In other words, the optical transition between the soliton 
and conduction bands rapidly develops as the soliton concentration increases. 
The energy positions of optical gaps of excitons are shown by the triangles at 



-10- 



the top of each figure. The lowest optical gap is about 0.7t, 0.9t, and l.Ot in 
Figs. 2(a), (b), and (c), respectively, and is slightly increasing. The almost 
constant behavior can be explained as follows. The HF order parameter Tj JiCr 
introduces an additional bond order, and thus increases the energy gap. In 
contrast, the lowest optical gap is a decreasing function of the soliton concen- 
tration in the free-electron case. Therefore, the decrease of the optical gap in 
the free electron case is suppressed by the increase of the energy gap in the 
presence of long-range Coulomb interactions. The optical gap of the transition 
between the continuum states is about 1.6t, l.7t, and 1.9t for the three figures. 
This quantity becomes larger with the concentration, due to the increase of 
number of states in the soliton band. 

We shall look at the optical spectra by HF and single-CI calculation in 
order to discuss exciton effects. They are shown in Fig. 3. The optical gap of 
the soliton-continuum exciton is 0.6t, 0.7t, and O.St in Figs. 3(a), (b), and (c), 
respectively. The optical gap of the intercontinuum exciton is lAt, lAt, and 
1.6t for the three figures. Here, we have regarded that the energy position of 
the lowest excitation where the component of the soliton-continuum exciton 
becomes smaller than 0.5, is the optical gap of the intercontinuum exciton. 
Both optical gaps decrease apparently from those of Fig. 2. This is due to the 
binding of an electron and a hole in the CI treatment. We also find that the 
soliton-continuum exciton has the larger total oscillator strength than that of 
Fig. 2. This is due to the one dimensionality, discussed in ref. [5]. 

There are many small structures in the optical spectra due to the finite 
system size. They could be removed by doing calculations for larger systems, 
as has been done in the paper [13] with reducing the dimension of the matrix 
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of CI excitations by means of the translational symmetry for the half-filled 
hamiltonian. But, the reduction of the system size is difficult for soliton lat- 
tice states owing to the periodicity of the system which changes with the soliton 
concentration. In order to pursue the change of optical excitation characters 
systematically for various combinations of parameters, we rather perform cal- 
culations for small system sizes around 100 and analyze numerical data as 
functions of the soliton concentration. In fact, when the optical gap and the 
ratio of the total oscillator strengths of the soliton-continuum exciton are plot- 
ted against the soliton concentration, (iV cl — N)/N, the plots are arrayed rather 
smoothly. We shall look at the data in Figs. 4 and 5. 

Figure 4 summarizes the optical gaps of the two kinds of excitons. They are 
calculated by the HF followed by the single-CI. Figures 4(a) and (b) are for U = 
2V = 2t and = At, respectively. The gaps of the intercontinuum exciton and 
soliton-continuum exciton are shown by closed and open squares, respectively. 
The optical gap of the soliton-continuum exciton is almost independent of the 
concentration, owing to the balance between the decreasing due to the soliton 
concentration change and the widening of the energy gap owing to the presence 
of long-range Coulomb interactions. The gap of the intercontinuum exciton 
increases rapidly when the concentration is larger than 2.5 percent. This is 
due to the fact that the number of states in the soliton band increases, and 
thus the energy gap between continuum states increases. It seems that these 
properties are common for the two Coulomb parameter sets. The optical gaps 
are larger for stronger Coulomb interactions. This is due to the larger bond 
order parameters which enhance the Peierls gap of the system. 

Figure 5 shows the ratio of the oscillator strengths of the soliton-continuum 
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exciton, plotted against the soliton concentration. The closed and open squares 
are the results of the HF-CI and HF calculations, respectively. The closed 
squares have the larger ratio than the open ones. This is one of exciton effects. 
When the concentration is near zero, the ratio varies almost linearly. This 
would be the natural consequence for low concentrations, because interactions 
among solitons are exponentially small and thus the ratio is proportional to 
the number of solitons. The increase near the zero concentration is slower for 
the stronger U and V of Fig. 5(b) than in Fig. 5(a). This would be due to 
that the soliton width is smaller for the stronger Coulomb repulsions and the 
portion of regions with nearly perfect dimerization strengths is larger. The 
increase of the ratio saturates at about 5 percent in Fig. 5(a) and at about 
7 percent in Fig. 5(b). The soliton-continuum exciton becomes like a free 
exciton at larger concentrations owing to the formation of the soliton band. 

4. Bipolaron lattice systems: Confinement ef- 
fects 

There exist only two kinds of degenerate conjugated polymers. They are trans- 
polyacetylene and pernigraniline base. All the other conjugated polymers have 
nondegenerate ground states. Therefore, it is also interesting to look at exciton 
effects in this type of polymers. The structures of nondegenerate conjugated 
polymers are generally complex, including aromatic rings and side chains. By 
this reason, the simple SSH-type models directly simulate only the structures of 
linear-chain polymers: trans- and cis-polyacetylenes. However, the SSH model 
with the non-zero BK term has the nondegenerate ground state which is one of 



the general properties described by the simple SSH-type models. Therefore, it 
is of general interests to look at excitonic effects on nondegenerate conjugated 
polymers. We shall consider the model Eq. (1) with 5 = 0.02. This small BK 
term gives rise to the large energy difference between the ground state and the 
excited state with the reversed phase of the bond alternation pattern, after 
calculations of adiabatic approximation with the full lattice relaxation. The 
Coulomb and lattice parameters are the same as in the previous section. 

We here note that there exist conjugated polymers, where phenylene rings 
are present, such as, polyaniline and poly(para-phenylenevinylene) (PPV). It 
is possible that a state with the reversed bond alternation phase is not present. 
However, we believe that only a small Sq makes a large energy difference be- 
tween the ground state and the state with the reversed bond alternation phase, 
and thus general properties of nonlinear excitations do not depend sensitively 
on whether there are states with much larger energies than that of the ground 
state. 

Figure 6 shows the static lattice configuration and the excess-electron den- 
sity distribution for N = 100, iV cl = 104, U = At, and V = It. Figure 6(a) 
displays the lattice configuration. The ground state with the positive bond 
variable is more stable than the state with the negative bond variable. There- 
fore, the region with the negative bond variable becomes smaller than in Fig. 
1(a), and two neighboring solitons come closer each other to form a bipolaron. 
There are two bipolarons in Fig. 6(a). The electron distribution pattern in 
Fig. 6(b) reflects the fact that two solitons are confined to form a bipolaron. If 
the confinement is more stronger, two peaks in the charge density distribution 
change into a single peak. 
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As we have done in the previous section, we shall calculate optical ab- 
sorption spectra by using the HF wavefunctions only as well as by performing 
single-CI calculations. Figure 7 shows the results of HF absorption (thin lines) 
and those of HF plus single-CI absorption (thick lines). The following param- 
eters are the same as in Fig. 6: U — 2V — 4i and iV = 100. The electron 
number is N e \ = 102 for Fig. 7(a) and N d = 104 for Fig. 7(b). The broadening 
7 = 0.05t is used. The major difference between the HF and HF-CI absorption 
is that the overall feature in the HF absorption shifts to lower energies in the 
HF-CI one. This is one of the exciton effects. 

A single bipolaron has two midgap states. Then, there are two midgap 
bands in the bipolaron lattice system. The number of excitons is three. We 
shall call them as follows: the photo-excited state from the upper bipolaron 
band to the conduction band as the "lower bipolaron-continuum exciton" , the 
exciton from the lower bipolaron band to the conduction band as the "upper 
bipolaron-continuum exciton" , and the exciton between the continuum states 
as the "intercontinuum exciton". In Fig. 7(a), the optical gaps of the lower 
bipolaron-continuum exciton, the upper bipolaron-continuum exciton, and the 
inter-continuum exciton are about 0.5t, 0.7t, and 1.6t, respectively. In Fig. 
7(b), they are about 0.6i, 0.9t, and l.Qt, respectively. The total oscillator 
strength of the upper bipolaron-continuum exciton is always much smaller 
than that of the other two excitons. This property is already seen in the 
free electron case and is independent of the magnitudes of U and V. As the 
concentration of the bipolarons increases, the oscillator strength of the lower 
bipolaron-continuum exciton enhances rapidly, and that of the intercontinuum 
exciton becomes smaller. 
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In order to analyze the concentration dependences systematically, we per- 
form calculations for several combinations of the system size N and the elec- 
tron number N c \. The electron number is always even, and then the bipolaron 
number is half of the excess-electron number N c \ — N. Figure 8 summarizes the 
optical gaps of three kinds of excitons (single-CI description), plotted against 
the excess-electron concentration: (N c \ — N)/N. We take two combinations 
of Coulomb parameters: (U, V) = (2t, It) [Fig. 8(a)] and (4t,2t) [Fig. 8(b)]. 
The optical gaps of the intercontinuum exciton, and upper (lower) bipolaron- 
continuum excitons are shown by the filled, and crossed (open) squares, respec- 
tively. The arrays of the plots behave smoothly, so it seems that the size effects 
are small. The optical gaps of the intercontinuum exciton and upper bipolaron- 
continuum exciton increase gradually as functions of the excess-electron con- 
centration. But, the increase of the lowest optical gap is suppressed, and 
it is nearly constant for concentrations larger than about 5 %. The similar 
behaviors have been seen in Fig. 3 for the soliton lattice case. 

Finally, we consider character of optical excitations by looking at the ra- 
tios of the total oscillator strengths of lower and upper bipolaron-continuum 
excitons. Figure 9 shows the results plotted against the excess-electron concen- 
tration. The squares are the data of the lower bipolaron-continuum excitons, 
and circles are the data of the upper ones. The open symbols are the data of 
the HF absorption, and the closed ones are for the HF-CI calculations. The 
ratios in the HF-CI absorption spectra are calculated as we have done in the 
soliton lattice case. Here, we assign the optical excitation to the exciton whose 
component is the largest among three excitons. The closed squares have the 
larger ratio than the open ones, due to exciton effects. The increase near the 
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zero concentration of the lower bipolaron-continuum exciton is suppressed for 
the stronger U and V of Fig. 9(b) than in Fig. 9(a). In other words, the 
increase of the ratio is steeper for weaker Coulomb interactions. This would 
be due to that the bipolaron width is smaller for the stronger Coulomb re- 
pulsions and the portion of regions with nearly perfect dimerization strengths 
is larger. The increase of the ratio nearly saturates at about 5 percent. The 
lower bipolaron-continuum exciton becomes like a free exciton at larger con- 
centrations. In contrast, the ratio of the upper bipolaron-continuum exciton is 
always small, and it is smaller than about 0.20. This smallness is already seen 
in the free electron case, and persists when there are Coulomb interactions. 
The same property in the optical spectra of the polaron and the bipolaron has 
been discussed in [14]. 



5. Summary 

We have looked at exciton effects on soliton and bipolaron lattice states in 
a model of the interacting electron lattice system with long-range Coulomb 
interactions. The Hartree-Fock approximation and the single-CI method have 
been used to obtain optical absorption spectra. We have discussed the follow- 
ing properties: 

(1) By comparison of the HF absorption with the HF-CI one, we have seen 
exciton effects which are similar to those, discussed for the half-filled systems 
in refs. [5,7]. The attraction between the excited electron and the remaining 
hole makes the excitation energy smaller when the correlations are taken into 
account by the single-CI. The oscillator strengths of the lower excited states 
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become relatively larger than in the HF calculations. 

(2) We have looked at variations of relative oscillator strengths of two or three 
kinds of excitons described by the single-CI. While the excess-electron concen- 
tration is small, the ratio of the oscillator strengths of the lowest exciton in- 
creases almost linearly. And, more than 80 percent of the oscillator strengths 
accumulate at the lowest exciton when the excess-electron concentration is 
larger than about 5 percent. It seems that this accumulation is very rapid as 
the concentration increases. 

6. Discussion 

In this paper, we have discussed general properties of excitons in degenerate 
and nondegenerate conducting polymers. We believe that the general proper- 
ties are well described by the present calculations, even when phenylene rings 
are present in polymer chains. Recently, there have been discussions on the 
optical responses [15] and nonlinear excitaions [16] in the PPV. The envelope 
of the polaron excitation has been experimentally determined [16] in the PPV 
system. Even though the optical absorption spectra of the PPV could be well 
explained by the noninteracting model such as the SSH-type Hamiltonian, we 
will discuss in the next paper [17] that exciton effects should be considered 
properly in order to explain several experimental features, such as, the enve- 
lope of the localized spin around the polaron [16], the magnitude of the optical 
gap [15], and the photoconductivity threshold [15]. We can also extend the 
present calculation to the polaron lattice state of the PPV by considering the 
phenylene ring structures. This is an interesting problem for further studies. 
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The variations of the optical gaps and the ratios of the oscillator strengths 
as functions of the concentration of the nonlinear excitations (solitons or bipo- 
larons) are rather smooth. It seems that how to systematically deal with such 
kinds of the numerical data obtained from numerical diagonalizations of finite 
systems is not fully investigated yet. We have looked variations as functions 
of the concentration, and have found the smooth variations. This fact was not 
stressed upon earlier, and plotting as a function of the concentration will be 
one of the useful methods for future investigations of systems with nonlinear 
excitations. 
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FIGURE CAPTIONS 



Fig. 1. (a) The dimerization order parameter, (— l) n (y n+ i — y n )/2, and (b) 
the excess-electron density, (p n -i + 2p„ + p n+ i)/4. The parameters are 5 = 0, 
U = At, V = 2t, N = 100, and N cl = 104. See the text for the other 
parameters. 

Fig. 2. The optical absorption spectra calculated with the HF wavefunctions 
for (a) (N,N el ) = (101, 102), (b) (100,102), and (c) (100,104). The broadening 
7 = 0.05t is used. The units of the abscissa are arbitrary. The component 
of the electron-hole excitation between the soliton and conduction bands are 
shown by dots which are connected with thin lines. The dots are plotted with 
the same abscissa as the absorption data. The minimum of the component is 
zero, and the maximum is unity. The energy positions of the optical gaps are 
shown by the triangles at the top of each figure. 

Fig. 3. The optical absorption spectra calculated with the HF plus single-CI 
wavefunctions for (a) (N,N cl ) = (101,102), (b) (100,102), and (c) (100,104). 
The broadening 7 = 0.05t is used. The units of the abscissa are arbitrary. The 
component of the electron-hole excitation between the soliton and conduction 
bands are shown by dots which are connected with thin lines. The dots are 
plotted with the same abscissa as the absorption data. The minimum of the 
component is zero, and the maximum is unity. The energy positions of the 
optical gaps are shown by the triangles at the top of each figure. 
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Fig. 4. The optical gaps in the single-CI of the "intercontinuum exciton" 
(filled squares) and of the "soliton-continuum exciton" (open squares), plot- 
ted against the soliton concentration. Coulomb interaction parameters are 
(U, V) = (2t, It) for (a), and (At, 2t) for (b). We note that two symbols at the 
5% concentration are the data for the combinations, (N, N c \) = (80, 84) and 
(120,126). 

Fig. 5. The ratio of the total oscillator strength of the "soliton-continuum 
exciton" as a function of the soliton concentration. The open squares are the 
data for the HF absorption, while the filled ones for the HF-CI absorption. 
Coulomb interaction parameters are (U,V) = (2t,lt) for (a), and (4t,2t) for 
(b). 

Fig. 6. (a) The dimerization order parameter, (— l) n (y n+ i — y n )/2, and (b) the 
excess-electron density, (p n _i + 2p n + p n+1 )/4. The parameters are 5 = 0.02, 
U = At, V = 2t, N = 100, and iV c i = 104. See the text for the other 
parameters. 

Fig. 7. The optical absorption spectra calculated with the HF (thin lines) 
and HF-CI (full lines) wavefunctions for (a) (N, N cl ) = (100, 102) and (b) 
(100,104). The broadening 7 = 0.05t is used. The units of the abscissa are 
arbitrary. 

Fig. 8. The optical gaps in the single-CI of the "intercontinuum exciton" (filled 
squares), of the "upper bipolaron-continuum exciton" (crossed squares), and of 
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the "lower bipolaron-continuum exciton" (open squares), plotted against the 
excess-electron concentration. Coulomb interaction parameters are (U, V) = 
(2t,l*) f o r (a), and (At, 2t) for (b). 

Fig. 9. The ratio of the total oscillator strength of the "upper (lower) 
bipolaron-continuum exciton" as functions of the excess-electron concentra- 
tion. The squares are for the lower exciton, and the circles are for the upper 
exciton. The open symbols are the data for the HF absorption, while the 
filled ones for the HF-CI absorption. Coulomb interaction parameters are 
(U, V) = (2t, It) for (a), and (At,2t) for (b). 
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